Parameter free scaling relation for nonequilibrium growth processes 
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We discuss a parameter free scaling relation that yields a complete data collapse for large classes 
of nonequilibrium growth processes. We illustrate the power of this new scaling relation through 
various growth models, as for example the competitive growth model RD/RDSR (random deposi- 
tion/random deposition with surface diffusion) and the RSOS (restricted solid-on-solid) model with 
different nearest-neighbor height differences, as well as through a new deposition model with tem- 
perature dependent diffusion. The new scaling relation is compared to the familiar Family- Vicsek 
relation and the limitations of the latter are highlighted. 
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The study of growing interfaces has been a very active 
field for many years pj|2j 0] . Many studies focus on the 
technologically relevant growth of thin films or nanos- 
tructures, but growing interfaces are also encountered in 
various other physical, chemical, or biological systems, 
ranging from bacterial growth to diffusion fronts. Over 
the years important insights into the behavior of nonequi- 
librium growth processes have been gained through the 
study of simple model systems that capture the most im- 
portant aspects of real experimental systems 0, [1| . 

In their seminal work, Edwards and Wilkinson inves- 
tigated surface growth phenomena generated by particle 
sedimentation under the influence of gravity [6(. They 
proposed to describe this process in (d + 1) dimensions 
by the following stochastic equation of motion for the 
surface height h(x, t), now called the Edwards- Wilkinson 
(EW) equation, 



dh(x,t) 
dt 



(1) 



where v is the diffusion constant (surface tension), 
whereas 77 (x, t) is a Gaussian white noise with zero mean 
and covariance (t?(x, t)r)(y, s)) = D5 d (x — y)S(t — s). 
Since Eq. [T] is linear\ it can be solved exactly by Fourier 
transformations 0, 0, 0|- Later, Family 0] discussed the 
random deposition (RD) and random deposition with 
surface relaxation (RDSR) processes. RD 0, 0] is one 
of the simplest surface growth processes. In this lat- 
tice model particles drop from randomly chosen sites 
over the surface and stick directly on the top of the 
selected surface site. Since there is no surface diffu- 
sion, the independently growing columns yield an un- 
corrected and never-saturated surface. The RDSR pro- 
cess is realized by adding surface diffusion which allows 
particles just deposited on the surface to jump to the 
neighboring site with lowest height. This diffusion step 
smoothes the surface and limits the maximum interface 
width W(t), defined at deposition time t as the stan- 
dard deviation of the surface height h from its mean 

value h: W(t) = {Jh — h) 2 ^ . Starting from an ini- 
tially flat surface, RDSR yields at very early times, with 



t < ti <~ 1 (we assume here that one layer is deposited 
per unit time), a surface growing in the same way as for 
the RD process since no (or only very few) diffusion steps 
occur in that regime. For t > t\ the width increases as a 
power law of time with a growth exponent before en- 
tering a saturation regime after a crossover time see 
Fig. Q] Both the saturation width W2 and the crossover 
time ti are powers of the substrate size L: 



Wo 



(2) 



where a is the roughness exponent and z is the dynami- 
cal exponent, with z = a/ [3. In his study Family noticed 
that the scaling exponents obtained through numerical 
simulations of the RDSR process agree with those ob- 
tained from the solution of the EW equation. The de- 
pendence of the growing interface on the substrate size 
L is summarized in the celebrated Family- Vicsek scaling 
relation [81 



W = L a f(t/L z ) 



(3) 



Combining this with the relations given in we see 
that this scaling relation mainly consists in shifting the 
crossover points for the different system sizes to a com- 
mon point with the new coordinates t' = t/L z ~ £/t 2 
and W = W/L a - W/W 2 . It is worth noting that the 
Family- Vicsek relation neglects the RD regime at early 
times and exclusively focuses on the two regimes con- 
nected by the crossover point at t = <2- 

The scaling ((2]) and ([3|) is generic for growing inter- 
faces and has been verified analytically, numerically, and 
experimentally in a large variety of systems. Various uni- 
versality classes have been identified which differ by the 
values of the scaling exponents. Thus the RDSR process 
belongs to the Edwards- Wilkinson universality class with 
the exponents (3—1/4 and z — 2 for a one-dimensional 
substrate. The RD process is in an universality class of 
its own which for a d=l substrate is characterized by the 
values (3=1/2 and z = 00. Other well known univer- 
sality classes, directly related to technologically relevant 
growth processes, are the Kardar-Parisi- Zhang (KPZ) Q 
and the conserved KPZ universality classes (lQ. [Tll|. 
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FIG. 1: Schematic plot of the interface width as a function of 
time for a typical deposition process. The early time behavior 
(with t < ti) is that of the random deposition (RD) process. 
For t > t\ correlated growth sets in. Finally, the finite system 
displays at late times t > t-2 a crossover to a saturation regime 
where the width of the interface remains constant, W = W2. 



In recent years the study of nonequilibrium growth sys- 
tems has shifted to more complex cases as for example 



competitive 



rowth models, 
23 



see e.g. 



H, El EI El Ei 



13, 0, El El El El, |23| . In a competitive growth model 
one considers a mixture of two different deposition pro- 
cesses where one of them takes place with probability p 
whereas the other takes place with probability 1 —p. One 
example is the RD/RDSR model [12| where the deposi- 
tion happens according to the RDSR rules with prob- 
ability p and to the RD rules with probability 1 — p. 
Whereas for p — 1 and p = only one of the processes 
is realized, for general values of p the mixture of the two 
processes leads to a crossover between the two regimes 
where the crossover time and width depend on the value 
of p (see Fig. [2^,). A similar dependence on system pa- 
rameters is also observed in the RSOS model [24| which 
exhibits a crossover from the RD regime to the important 
KPZ universality class. In this model new particles are 
incorporated into the growing surface only if the height 
differences between the deposition site and its neighbor- 
ing sites remain smaller than some maximum height S. 
As discussed in (2f| and shown in Fig. [3^, the crossover 
time and width depend on the value of S. 

In simple growth processes the random deposition 
regime is restricted to very early times. This is fun- 
damentally different in more complex systems where 
the initial regime can extend over very l arg e times 
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ready mentioned, the Family- Vicsek scaling relation ^ 
assigns a new set of coordinates to the second crossover 



point. This does however not yield a complete data col- 
lapse for growth processes with two crossover points if 
one considers systems of different sizes. For the com- 
petitive growth models some phenomenological scaling 
relations have been proposed in the past, but these mod- 
ified scaling relations also only allow a partial collapse of 
the different curves [H, El El HI . 

However, a scaling relation leading to a complete data 
collapse of all curves obtained for different system sizes 
and different values of the system parameters can indeed 
be obtained for any growth system that exhibits two dif- 
ferent crossover points. This data collapse is achieved in 
a two-step process. First we translate all curves in a log- 
log plot such that the first crossover point is now located 

at the origin. This is achieved by plotting log (jfcj as a 
function of log ( 4- ) • In the second step we rescale both 



axes by the common scale factor A = 1 / log ^ j such 

that in the log-log plot the second crossover point is fixed 
at the rescaled width logW^ = 1- This isotropic reseat- 
ing, which conserves the slope of the region between the 
two crossover points, makes that the length of the line 
connecting the two crossover points is the same for all 
curves, and a complete data collapse, encompassing all 
three regimes, follows. Our proposed scaling relation can 
be cast in the following equation: 



F 



log (fr 



(4) 



where F(x) is a scaling function. Introducing A 
1/ log [y^J , we can rewrite this as 



W 



(5) 



with a new scaling function G(y). As shown in Fig. [2] 
and [3] for the RD/RDSR and RSOS processes, the pro- 
posed scaling relation yields a complete data collapse for 
different system sizes L and different values of the sys- 
tem parameters. This perfect scaling behavior should be 
compared with the incomplete scaling proposed in the 
literature 0, El El HI- 

Obviously, the scaling relation (O is of universal use 
in growth systems with two crossover points and replaces 
the Family- Vicsek relation in these systems. This class 
of systems encompasses competitive growth models, but 
also the simple growth systems, for which the Family- 
Vicsek relation has been proposed originally, belong to 
this class. It is also worth noting that the properties of 
the different models only enter in our relation ([5]) im- 
plicitly through the dependence of the positions of the 
crossover points on the different system parameters. 

We can also state the new scaling relation in an al- 
ternative way which makes the difference to the Family- 
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FIG. 2: (a) Log-log plot of the surface width vs time for 
the RD/RDSR process in systems with different sizes L and 
different probabilities p. (b) A complete data collapse of all 
data sets is achieved when using the scaling relation ((5j). The 
curves shown in the left panel completely fall one on top of 
the other and are no longer distinguishable. The axis labels 
are t' = t x /t x and W' = W x /W x , see main text. The dashed 
lines indicate the expected slopes in the RD and EW regimes. 
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FIG. 3: The same as Fig. H but now for the RSOS model 
with different values S of the maximal nearest-neighbor height 
difference. Again a perfect data collapse is achieved when 
using the scaling relation ((5j). The dashed lines indicate the 
expected slopes in the RD and KPZ regimes. 



W x = W$ G 



(7) 



where F and G are again scaling functions. This scaling 
relation is completely equivalent to the first one, only 
the scales are shifted. In fact, the relation |7]) allows 
a direct comparison with the Family- Vicsek relation |[3| : 
recalling that the scaling behaviors of W% and t 2 are given 
by the relations @, we immediately see that we recover 
the Family- Vicsek relation by setting A = 1. This nicely 
shows that it is the isotropic rescaling in the log-log plot 



by the factor 



that ultimately is responsible for 



the success of the new scaling relation. 

Competitive growth models have the peculiar feature 
that at every deposition one has to decide which of the 
two deposition rules is followed by the newly added par- 
ticle. We propose in the following a deposition model 
with similar features as the competitive growth models, 
but where the competition is intrinsic and governed by 
the value of the substrate temperature. This is a much 
more realistic situation, especially since in the growth of 
thin films and nanostructures the substrate temperature 
is an important parameter that shapes to a large extend 
the morphology of growing structures [2r3 |. 

The deposition model discussed in the following is 
based on Family's original RDSR process and differs 
from this model by the diffusion step. In the RDSR pro- 
cess a particle deposited on the surface is allowed to jump 
to one of the neighboring sites if this site has a lower 
height than the site of deposition. In our model we as- 
sign an energy -E(x, t) — f h(x, t) to the column x where 
h(x,t) is the height of that column at time t. The con- 
stant / can be thought to be the gravitation constant. 
Starting from an initially flat substrate, particles are de- 
posited on randomly chosen sites and then allowed to 
diffuse locally after deposition. For a diffusion step tak- 
ing place at time t, we select one of the neighboring sites 
y at random and accept the jump with the temperature 
and time dependent (Metropolis like) probability 



1 , if E(y,t) <£(x,t) 

e -[E(y,t)-E(x,t)]/k B T 

= e -f[h(y,t)-h(x,t)]/k B T otherwise. 



Vicsek relation more transparent. Indeed, a complete 
collapse can also be achieved when first moving the sec- 
ond crossover point to the origin in a log- log plot, yielding 
the relation 



F 



^ (ft 



(6) 



In the following we choose units thus that f/ks = 1 
where ks is the Boltzmann constant. 

In contrast to the original RDSR model, we have in 
the present model a non- vanishing probability that a de- 
posited particle jumps to a neighboring site with a higher 
height than the deposition site. We assume this jump to 
be thermally activated and to depend on the tempera- 
ture T of the substrate. As we discuss elsewhere [13] > 
the substrate temperature is a parameter that allows the 
study of novel properties of growing interfaces. 
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In Fig. [4ji we show the temporal evolution of the width 
for various temperatures and system sizes. As for the 
RDSR process one distinguishes three regimes separated 
by two crossover points: a RD regime, followed by a EW 
regime, with a final crossover to the saturation regime. 
In contrast to the RDSR process, the random deposition 
process is not confined to the very early time regime t < 1 
but can extend to larger times. In fact, the crossover time 
t\ between the random deposition and the EW regimes is 
shifted to higher values for increasing temperatures and 
diverges in the limit of infinite temperatures. Of special 
interest is that the surface widths shown in Fig. [4^ can be 
directly obtained from the Edwards- Wilkinson equation 
([1]) for a system of size L with a temperature dependent 
diffusion constant. For example, for the temperatures 
shown in Fig.gh., the values of v are: v(T — 100) = 0.005, 
u(T = 10) = 0.04, and u(T = 1) = 0.18. 
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FIG. 4: (a) Log-log plot of the surface width vs time for the 
temperature dependent deposition model. Systems of differ- 
ent sizes at different temperatures are shown, (b) A complete 
data collapse of all data sets is also achieved for this model. 
The dashed lines indicate the expected slopes in the RD and 
EW regimes. 

We first check in Fig.QJ) that also for the present model 
the scaling relation ((SJ yields the full data collapse. Due 
to the simplicity of the model, we can obtain the full 
information on the location of the two crossover points 
27| . In this way we find that, as usual, only the crossover 
to the saturation regime depends on the system size. In 
addition, the coordinates of both crossover points display 
a linear dependence on the substrate temperature. Tak- 
ing these observations into account, we can rewrite the 
scaling relation {5} in the form 



W x = (a + bT) X/2 G 



(a + bry 



(8) 



of the EW universality class. Eq. ([8j) directly reveals 
for our model the dependence of the generalized scaling 
relation on the system size and on the temperature. 

In conclusion, we have presented in this letter a pa- 
rameter free scaling relation that yields a complete data 
collapse for large classes of nonequilibrium growth pro- 
cesses with two crossover points. Examples include all 
simple growth processes as well as more complex grow- 
ing interfaces as encountered for example in competitive 
growth systems. A deposition model with temperature 
dependent diffusion allows us to discuss the dependence 
of the scaling relation on the relevant system parameters. 



where a = 0.59 and b = 0.29, with A(L) = 1 /log (c L a ) 
and c = 0.55, whereas a = 1/2 is the roughnes exponent 
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